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基于 神经 网 络 的 准 实时 单 站 电离 层 TEC 反 演 


卢 伟 俊 12， 马 冠 一 
(1. 中 国 科 学 院 国 家 天 文 台 ， 北 京 100101; 2. 中 国 科 学 院 大 学 ， 北 京 100049) 

摘要 : 使 用 全 球 导航 卫星 系统 (GNSS) 进行 电离 层 总 电子 密度 (TEC) 反 演 时 ， 单 站 相对 于 多 站 观测 网 是 
一 种 灵活 简便 的 方法 。 基 于 人 工 神经 网 络 ， 本 文 提出 了 一 种 准 实时 单 站 电离 层 TEC 反 演 的 方法 。 在 这 种 方法 中 ， 
上 一 个 时 段 的 硬件 偏差 被 作为 初 值 并 随 着 观测 值 调整 ， 同 时 电离 层 TEC 也 被 准 实时 地 反 演 。 为 了 对 这 种 方法 进 
行 详细 的 评估 ， 通 过 位 于 中 国 的 单 站 ， 四 天 磁 静 日 的 TEC 被 分 别 采 用 提出 的 方法 与 经 典 的 最 小 二 乘 球 谐 函数 法 
反 演 ， 其 中 硬件 偏差 和 TEC 的 参考 值 通过 附近 的 多 站 网 得 到 。 在 另 一 个 测试 中 ， 通 过 位 于 欧洲 的 单 站 ， 一 次 电 
离 层 暴 事件 及 前 后 平静 日 的 TEC 也 被 分 别 通 过 上 述 方法 反 演 。 在 磁 静 日 ， 估 计 的 硬件 偏差 在 整体 上 相对 于 最 小 
二 乘 球 谐 函 数 法 更 接近 于 参考 值 ， 反 演 的 TEC 也 更 接近 于 参考 值 。 电 离 层 暴 时 ， 两 种 方法 反 演 的 TEC 也 具有 很 
好 的 一 致 性 ， 神 经 网 络 法 估计 的 硬件 偏差 与 磁 静 日 的 硬件 偏差 更 接近 。 结 果 表 明 ， 提 出 的 神经 网 络 法 相 比 最 小 二 
乘 球 谐 函数 法 具有 较 高 的 精度 
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电离 层 是 位 于 地 面 上 约 60-1000 和 干 米 高 度 的 大 气 区 域 , 其 中 部 分 中 性 粒子 被 电离 为 自由 电子 和 离子 由。 
因此 ， 电 离 层 是 一 种 色散 介质 并 引起 无 线 电 信号 的 群 延迟 和 相位 超前 ， 影 响 程度 与 信号 频率 及 电离 层 电子 
密度 有 关 。 为 了 评估 电离 层 对 无 线 电 信和 号 传播 的 时 延 影响 ， 一 些 描述 电离 层 电子 密度 的 理论 模型 或 经 验 模 
型 被 提出 DB]。 理 论 模型 虽然 能 反映 电离 层 变 化 的 物理 机 制 ， 但 在 精度 上 还 无 法 满足 电离 层 延迟 改正 的 要 
求 。 而 经 验 模型 和 真实 电离 层 之 间 也 存在 着 无 法 忽视 的 偏差 ， 只 能 改正 部 分 电离 层 延迟 。 因 此 ， 通 过 真实 
观测 对 电离 层 进行 建 模 已 经 成 为 一 种 必要 的 手段 。 一 些 观 测 仪器 例如 测 高 仪 和 非 相干 散射 雷达 已 经 被 用 于 
【em 探测 电离 层 的 电子 密度 ， 然 而 这 些 仪器 只 能 得 到 测 站 上 空 的 电离 层 剖 面 而 无 法 进行 大 规模 的 电离 层 观 测 。 
© 随 着 全 球 导 航 卫星 系统 (Global Navigation Satellite System, GNSS) 的 运行 ， 电 离 层 观测 可 以 通过 GNSS 

多 频 无 线 电信 号 进行 。 由 于 GNSS 是 基于 倾斜 轨道 卫星 星座 连续 运行 的 ， 不 同 的 时 空 尺 度 的 电离 层 变化 都 

可 以 通过 GNSS 进行 探测 。 

利用 GNSS， 电 离 层 电子 密度 的 积分 测量 值 即 总 电子 含量 (Total Electron Content, TEC) 被 作为 参数 

指标 进行 电离 层 反 演 。 不 同 频 率 无 线 电信 号 通过 卫星 和 接收 机 硬件 通道 时 产生 的 硬件 偏差 作为 系统 误差 必 

须 被 剔除 ，GNSS 多 站 观测 网 通常 被 设计 用 于 电离 层 硬 件 偏 差 的 估算 和 电离 层 参 数 的 反 演 。 尽 管 如 此 ， 单 

站 电离 层 反 演 作为 Gi A i 究 价 值 。 单 站 电离 层 反 演 的 算法 主要 有 网 格 

法 、 卡 尔 曼 滤波 法 和 多 项 式 法 四。 多 项 式 法 中 的 多 项 式 可 以 是 二 次 多 项 式 、 三 角 级 数 或 球 谐 函 数 钻 。 一 般 

来 说 二 次 多 项 式 和 三 角 级 数 只 适合 较 小 范围 的 反 演 ， 球 谐 函 数 通 过 改变 阶 数 和 次 数 能 更 好 地 描述 电离 层 在 
不 同 空间 尺度 的 变化 。 

对 于 电离 层 反 演 中 的 拟 合 方法 ， 最 小 二 乘 和 卡尔 曼 滤 波 通常 会 被 采用 。 随 着 人 工 神经 网 络 〈Artificial 
Neural Network, ANN) 被 提出 ， 它 很 快 就 被 用 于 进行 电离 层 反 演 器。 一 个 基于 误差 反馈 的 前 向 传播 神经 网 
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络 能 够 以 任意 精度 逼近 一 个 连续 的 实 函 数 。 这 种 拟 合 是 通过 迭代 的 方式 进行 的 ， 各 参数 对 应 的 权 值 可 以 动 
态 地 调整 。 因 此 这 种 动态 拟 合适 合 于 电离 层 的 准 实 时 反 演 。 

基于 人 工 神经 网 络 ， 我 们 提出 一 种 准 实 时 的 单 站 电离 层 反 演 方法 。 首 先 对 电离 层 反 演 原理 和 提出 的 方 
法 进行 了 推导 。 接 下 来 利用 位 于 北纬 30.28° 东 经 109.46? 的 单 站 在 2014 年 春分 、 夏 至 、 秋 分 和 冬至 的 GNSS 
数据 进行 了 实验 ， 其 中 位 于 东经 110°? 子 午 线 的 18 个 站 被 作为 验证 站 。 在 另 一 个 实验 中 ，2016 年 10 月 13 


日 电离 层 烘 事件 的 数据 也 被 用 于 电离 层 反 演 ， 并 且 结 果 与 前 后 磁 静 日 的 结果 进行 了 比较 。 最 后 一 部 分 是 对 
文章 的 总 结 


1 电离 层 反 演算 法 


1.1 球 谐 函 数 拟 合 
GNSS 卫星 以 两 个 不 同 的 工 波段 频率 播发 信号 , 使 用 双 频 接收 机 , 伪 距 和 载波 相位 观测 值 可 以 被 接收 : 
P;=p+t+c(6,-—6,)-I;+T +e (1) 
Li =p +c(ô;— ôs) +l; +T +N; +E (2) 


fi 

其 中 i 为 1 或 2 分 别 对 应 所 或 f2 频 率 ，P; 和 Li 分 别 对 应 f 频率 的 伪 距 和 载波 相位 观测 值 ，p 是 卫星 和 接收 机 

间 的 真实 距离 ，c 表示 光速 ，654 和 6, 分 别 是 接收 机 和 卫星 钟 差 ， 并 对 于 伪 距 包含 群 延迟 ， 对 于 载波 相位 包 

含 相 位 超前 ;五 是 对 应 广 频 率 的 电离 层 延迟 ，7 是 对 流 层 延迟 ，Ni; 为 L; 的 整 周 模糊 度 ，e 表 示 观 测 值 的 残 差 
项 。 

由 不 同 频 率 观 测 值 的 差分 ， 沿 着 传播 路 径 的 斜 向 TEC (Slant TEC, STEC) 可 以 分 别 通过 伪 距 和 载波 

相位 观测 值 求 出 : 


Fafa)? 
SRC se Pa) G3) 


(fif2)? cla cla , cN1 cN2 
STEC, = pea Se et Op ek eae 4 
E 40.3(f2-f3) 万 f2 fi f2 ) (4) 


考虑 到 伪 距 包含 较 大 观测 噪声 ， 且 载波 相位 含有 未 知 的 整 周 模糊 度 ， 高 精度 STEC 可 以 通过 载波 相位 
平滑 伪 距 求 得 吕 ; 
STEC(t) = w,STECp(t) + (1 — w,) [STEC(t — 1) + STEC,(t) — STEC,(t—1)] (t > 1) (5) 
其 中 t 代 表 历 元 ，w, 表 示 权 重 并 可 以 与 卫星 高 度 角 有 关 。 对 于 第 一 个 历 元 ，STEC 与 STECp 相 等 。 
电离 层 反 演 即使 用 STEC 得 到 反 演 区 域内 沿 重 直方 向 的 TEC。 通常 假设 电离 层 的 所 有 电子 集中 于 一 个 
无 限 薄 的 单 层 上 ， 为 此 引入 了 如 图 1 所 示 的 电离 层 单 层 模型 中 。 其 中 点 S$、R 和 0 分 别 表示 卫星 、 接 收 机 
和 地 心 。 卫 星 和 接收 机 视线 与 电离 层 单 层 的 交点 被 定义 为 穿刺 点 Clonospheric Piercing Point, IPP), HZ 


纬度 等 于 对 应 的 星 下 点 〈Sub-Ionospheric Point, SIP) 的 经 纬度 。 因 此 TEC 可 以 根据 图 1 所 示 的 几何 关系 
得 出 : 

TEC = (STEC — bs—b;)'cosx (6) 

X 三 arcsin( E) (7) 


FLA b Mb AIE PEMA AS A mE, A E 分 别 是 天 顶 角 和 高 度 角 ，R。 是 地 球 的 平均 半径 ， 互 是 
电离 层 单 层 的 高 度 。 
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图 1 电离 层 单 层 模型 


Fig. 1 Ionospheric single layer model 


在 另 一 方面 ，TEC 在 建立 的 模型 中 可 以 被 认为 是 一 个 球 谐 函数 ， 参 数 为 IPP 对 应 的 余 纬 和 日 地 坐标 系 


T FRAR, 

- TEC = a Er o (Anm cos MP + Bum sinmg )Pnm( cos 0) (8) 
=) 其 中 g 和 gq 分 别 为 余 纬 和 日 地 坐标 系 下 的 经 度 ， 对 于 单 站 ， 通 常 取 阶 数 和 次 数 为 5 CN = M = 5);， Pa IE 
O ”一 化 的 勒 让 德 函 数 ，A4nm 和 Bm 为 球 谐 函 数 的 系数 ， 并 可 以 通过 最 小 二 乘 或 其 它 拟 合 方法 求 得 。 

1O ”1.2 人 工 神经 网 络 

© 人 工 神经 网 络 由 许多 神经 元 组 成 ， 这 些 神经 元 组 成 了 一 个 输入 层 、 一 个 输出 层 和 若干 隐藏 层 。 通 常 一 
A 个 隐藏 层 是 足够 的 。 一 个 基于 误差 反馈 的 多 层 前 向 传播 神经 网 络 能 够 以 任意 的 精度 趋 近 于 一 个 连续 实 函数 
名 (Slo 如 图 2 所 示 为 用 于 准 实时 电离 层 反 演 的 神经 网 络 的 示意 图 ， 其 中 大 表示 卫星 和 接收 机 的 数目 。 由 于 在 


> 电离 层 反 演 时 ， 卫 星 和 接收 机 的 硬件 偏差 需要 同时 被 估计 并 被 剔除 ， 额 外 的 网 络 被 设计 用 于 硬件 偏差 的 佑 
计 。 输 入 参数 为 STEC 对 应 的 余 纬 、 日 地 坐标 系 下 的 经 度 以 及 高 度 角 。 输 出 参数 为 式 (8) 所 示 的 球 谐 函 
数 的 系数 以 及 硬件 偏差 , 那么 观测 站 上 空 的 TEC 可 以 通过 球 谐 函 数 的 系数 求 得 。 输 入 层 的 神经 元 除 硬 件 偏 
差 的 初 值 外 ， 还 有 IPP 的 余 纬 % 和 日 地 坐标 系 下 的 经 度 9 以 及 天 项 角 X。 和 输出 层 包 括 硬件 偏差 以 及 球 谐 函数 
TT 系数 。 隐 藏 层 的 节点 数 是 一 个 小 于 200 的 随机 值 。 输 出 层 的 值 通过 动态 调整 的 权 值 进行 调整 。 权 值 通过 梯 
度 下 降 算法 根据 观测 值 与 模型 估计 值 的 残 差 求 出 叫 ; 

E = 1/2- [VTEC/cos (x) — (STEC 一 DCBD)]? (9) 
根据 式 9)， 以 每 15 分 钟 时 段 内 的 各 卫星 对 应 的 STEC 为 训练 数据 集 ， 在 对 神经 网 络 进行 训练 的 过 
不 断 地 调整 权 值 和 输出 层 节 点 ， 最 终 使 式 (9) 收敛 。 由 于 电离 层 的 变化 具有 连续 性 ， 下 一 个 时 段 可 
用 当前 时 段 训 练 的 权 值 和 输出 层 节 点 作为 初始 值 , 并 继续 用 新 的 数据 集 进 行 训练 , 以 达到 快速 收敛 的 目的 。 
对 于 第 一 个 时 段 ， 前 一 天 的 硬件 偏差 被 作为 初 值 ， 之 后 每 一 个 时 段 都 用 前 一 个 时 段 的 硬件 偏差 作为 初 值 。 
输入 值 经 过 非 线性 变换 进入 输入 层 再 经 过 隐 单 元 逐 层 处 理 , 传 向 输出 层 , 输出 层 判断 估 值 与 观测 值 的 残 差 ， 
大 于 阔 值 则 转 入 反 向 传播 ， 通 过 梯度 下 降 算 法 修改 各 神经 元 的 权 值 ， 直 至 输入 满足 条 件 。 由 于 不 同 卫星 对 
应 的 观测 值 包含 非 均 一 的 观测 误差 LOU， 而且 电离 层 受到 的 扰动 也 会 随 着 时 间 动 态 变 化 ， 因 此 神经 网 络 是 
一 种 良好 的 准 实时 估计 算法 。 

我 们 基于 图 2 所 示 的 神经 网 络 架 构 使 用 MATLAB 语言 编写 了 准 实 时 单 站 电离 层 反 演 程序 。 该 程序 每 
15 分 钟 通过 实时 接收 到 的 GNSS 数据 进行 一 次 电离 层 反 演 ， 得 到 一 组 球 谐 函数 系数 ， 并 更 新 硬件 偏差 。 使 
用 3.4 Ghz 的 CPU 和 15 GB 的 RAM 的 工作 站 测试 后 证 实 每 一 次 电离 层 反 演 可 以 在 一 分 钟 左 右 完成 ,因此 
可 以 实现 准 实时 的 单 站 电离 层 反 演 。 
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图 2 人 工 神 经 网 络 示意 图 


Fig. 2 Schematic diagram of artificial neural network 


2 磁 静 日 数据 验证 


2.1 测试 观测 网 和 数据 处 理 方 法 

如 图 3 所 示 为 GNSS 观测 站 的 分 布 , 其 中 位 于 北纬 30.28° 东 经 109.46* 的 观测 站 为 测试 站 ,其 余 位 于 东 
经 110° 子 午 线 附近 的 18 个 站 为 验证 站 。 为 了 对 提出 的 准 实时 神经 网 络 法 RTANN) 进行 详细 的 评估 ， 我 
们 使 用 测试 站 的 全 球 定位 系统 (Global Positioning System, GPS) 的 双 频 观测 数据 分 别 用 准 实 时 神经 网 络 
法 和 最 小 二 乘 球 谱 函数 法 (LSM) 反 演 电离 层 TEC。 同 时 ， 使 用 18 个 验证 站 的 GPS 双 频 观测 数据 采用 1° 
x le 的 自 适 应 网 格 法 估计 硬件 偏差 05， 再 使 用 2。x 2 的 网 格 法 每 15 分 钟 得 到 电离 层 TEC 03， 并 将 它们 作 
为 参考 值 与 准 实时 神经 网 络 法 和 最 小 二 乘 球 谐 函 数 法 得 到 的 结果 进行 比较 。 自 适应 网 格 法 已 经 在 文献 [12] 
中 应 用 于 电离 层 暴 时 的 硬件 偏差 估计 ， 因 此 本 文 使 用 该 方法 得 到 硬件 偏差 参考 值 。 除 此 之 外 ， 尽 管 多 站 网 
的 反 演 结果 也 含有 误差 ， 但 在 观测 站 分 布 合理 的 情况 下 ， 多 站 网 的 反 演 结果 比 单 站 有 更 高 的 精度 和 更 好 的 
稳定 性 ， 因 此 作为 参考 值 。 在 实验 中 ， 我 们 选择 了 2014 年 位 于 春分 、 夏 至 、 秋 分 和 冬至 观测 数据 ， 这 4 
天 均 为 磁 静 日 。 历 元 间隔 为 30 秒 ， 截 止 高 度 角 选取 20"， 电 离 层 单 层 高 度 设 为 400 km. 
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图 3 110°E 子午 线 的 GNSS 观测 站 的 分 布 
Fig. 3 Distributions of the GNSS observation stations along 100°E meridian 


2.2 估计 的 硬件 偏差 
由 于 在 电离 层 反 演 中 ，] 


a 


ChinaxX iv 
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星 和 接收 机 的 硬件 偏差 是 作为 整体 耦合 在 一 起 的 ， 我 们 采用 零 均 值 处 理 分 离 


卫星 和 接收 机 硬件 偏差 ， 即 对 每 一 个 卫星 硬件 偏差 减 去 它们 的 均值 并 将 均值 并 入 接收 机 硬件 偏差 。 图 4 为 
通过 准 实时 神经 网 络 法 和 最 小 二 乘 球 谐 函数 法 估计 的 卫星 硬件 偏差 与 参考 值 的 差 值 ，(a)-(d) 分 别 为 春分 、 
夏至 、 秋 分 和 冬至 的 结果 。 在 这 四 天 中 ， 通 过 准 实时 神经 网 络 法 得 到 的 卫星 硬件 偏差 与 参考 值 的 差 值 在 整 
体 上 比 最 小 二 乘 球 谐 函数 法 更 小 。 提 升 的 原因 在 于 神经 网 络 的 引入 使 得 硬件 偏差 可 以 动态 调整 。 这 样 减 小 


了 观测 噪声 以 及 电离 层 扰动 对 反 演 结果 的 影响 。 另 外 有 极 少数 卫星 的 硬 人 
由 零 均 值 处 理 造 成 的 。 还 有 部 分 卫星 在 测试 时 处 于 不 可 用 状态 ， 因 此 其 结果 在 图 中 没有 显示 。 经 过 计算 ， 


通过 准 实时 神经 网 络 法 估计 的 硬件 偏差 在 四 天 的 均 方 根 误差 分 别 为 0.3 
球 谐 函数 法 的 结果 为 0.84、0.61、0.67、0.68 ns。 可 以 发 现 春 分 时 的 提 刊 


F 偏 差 与 参考 值 的 偏差 较 大 ， 这 是 


8、0.25、0.37、0.31 nso | 


j 最 小 二 乘 


| 最 为 明显 ， 这 与 春分 时 


E BJZ TEC 


值 较 大 ， 对 硬件 偏差 估计 的 噪声 干扰 较 大 有 关 ， 这 将 在 下 文 进一步 分 析 。 
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图 4 不 同方 法 估计 的 卫星 硬件 偏差 与 参考 值 的 差 值 


Fig. 4 Deviations between instrumental biases of satellites estimated by different methods and the references 


2.3 反 演 的 TEC 
图 $ 为 通过 准 实时 神经 网 络 法 和 最 小 二 乘 球 谐 函 数 法 反 演 的 北纬 30.28" 东 经 109.46" 单 站 位 置 的 电离 层 
TEC，(a)-(d) 分 别 为 春分 、 夏 至 、 秋 分 和 冬至 的 结果 。 参 考 值 是 通过 18 个 站 采用 网 格 法 得 到 的 ， 并 且 每 15 
分 钟 独 立 反 演 。 可 以 发 现 图 5(a) 中 春分 时 的 TEC 明显 大 于 另外 三 天 ， 因 此 通过 式 (6) 与 TEC 同时 估计 的 
硬件 偏差 受 噪声 影响 也 更 大 。 将 单 站 反 演 的 TEC 与 参考 值 对 比 后 可 以 发 现 , 通过 准 实时 神经 网 络 法 反 演 的 
TEC 比 最 小 二 乘 球 谐 函数 法 更 接近 于 参考 值 ， 尤 其 是 正午 和 午夜 。 正 午 提升 较 明 显 是 因为 正午 时 TEC 值 
较 大 ， 而 神经 网 络 的 加 入 使 得 反 演 受 噪声 影响 更 小 。 而 午夜 是 因为 此 时 TEC 值 较 小 ， 最 小 二 乘 球 谐 函数 法 
的 反 演 结果 受 硬件 偏差 的 影响 较 大 ， 而 准 实时 神经 网 络 法 估计 的 硬件 偏差 更 准确 。 另 外 ， 准 实时 神经 网 络 
反 演 的 TEC 也 不 平滑 , 我 们 推测 这 可 能 与 电离 层 扰 动 有 关 ， 因为 大 部 分 波动 趋势 都 与 参考 值 相 对 应 。 Ait 
算 , 准 实时 神经 网 络 法 的 均 方 根 误差 分 别 为 2.76、1.00、0.99、1.73 TECU, 最 小 二 乘 球 谐 函 数 法 的 为 3.46、 
2.32、0.94、1.46 TECU。 因 此 ， 通 过 实时 神经 网 络 法 反 演 的 TEC 比 最 小 二 乘 球 谐 函 数 法 更 准确 。 
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Fig. 5 TEC at 30.28°N and 109.46°E derived by different methods 
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我 们 分 别 通 过 准 实时 神经 网 络 法 和 最 小 二 乘 球 谐 函 数 法 通过 单 站 反 演 了 2016 年 10 A 13 日 电离 层 暴 
的 TEC。 如 图 6 所 示 为 GNSS 观测 站 的 分 布 ， 其 中 位 于 北纬 40.52° 西 经 3.09° 的 观测 站 为 测试 站 ， 其 余 位 
F OFFRA 40°N 纬 线 附 近 的 7 个 站 为 验证 站 。 根 据 文献 [14]， 暴 的 起 始 发 生 于 12 日 2100 UT， 初 相 从 
2300 UT 持续 至 13 日 的 0600 UT, 接 下 来 主 相 持续 了 18 小 时 ; 并 且 使 用 图 6 的 7 个 验证 站 观测 到 13 日 白 
天 的 正 暴 及 夜间 的 负 暴 ， 结 合 测 高 仪 观测 到 的 电离 层 抬 升 证 实 了 引起 这 个 电离 层 暴 的 主要 物理 机 制 是 扰动 
中 性 风 。 
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图 6 0 子午 线 和 40°N 纬 线 的 GNSS 观测 站 的 分 布 
Fig. 6 Distributions of the GNSS observation stations along 0° meridian and 40°N parallel 


如 图 7 所 示 为 不 同方 法 估计 的 2016 年 10 月 8 日 至 17 日 的 卫星 硬件 偏差 相对 于 10 天 日 均值 的 均 方 根 
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F 偏 差 相 对 于 磁 静 日 变化 较 大 。 但 实际 上 硬件 偏差 取决 


于 仪器 本 身 以 及 所 处 的 外 部 环境 , 只 是 由 于 反 演 过 程 中 与 TEC 耦合 在 一 起 才 受 到 了 电离 层 变 化 的 影响 。 而 
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图 7 不 同方 法 估计 的 卫星 硬件 偏差 相对 于 日 均值 的 均 方 根 误差 


准 实时 神经 网 络 法 的 均 方 根 误差 在 这 10 天 中 均 小 于 最 小 二 乘 球 谐 函数 法 ， 特 别 是 在 13 日 暴发 生 期 间 。 这 
说 明 准 实时 神经 网 络 法 估计 的 硬件 偏差 受 电离 层 变 化 的 影响 更 小 ， 因 此 优 于 最 小 二 乘 球 谐 函 数 法 。 


Fig. 7 Root mean square error of instrumental biases of satellites estimated by different methods and the daily averages 


如 图 8 所 示 为 不 同方 法 反 演 的 2016 年 10 月 8 


日 至 17 日 


的 TEC。 相 比 于 8-11 日 磁 静 日 的 TEC，13 


日 白天 的 TEC 表现 为 明显 的 正 暴 效应 ，TEC 增 大 率 超 过 了 100%， 而 从 13 日 夜间 开始 表现 为 明显 的 负 暴 


效应 ， 并 持续 至 15 日 。 另 外 ， 反 演 的 TEC 与 文献 [14] 的 结果 相符 ， 并 且 与 赤道 向 扰动 中 性 风 的 物理 机 制 


相符 。 暴 时 的 赤道 向 中 性 风 可 以 抬升 电离 层 ， 同 时 也 会 引起 电离 层 的 成 分 变化 。 白 天 由 于 电离 层 较 低 高 度 


光电 离 作 用 的 补充 使 得 正 暴 发 生 ， 而 夜间 光电 离 作 用 消失 同时 离子 结合 率 上 升 使 得 负 暴 发 生 。 另 外 ， 比 


较 准 实时 神经 网 络 法 和 最 小 二 乘 球 谐 函数 法 的 反 演 结果 ， 发 现 两 种 方法 反 演 的 TEC 也 具有 很 好 的 一 致 性 。 
另外 ， 准 实时 神经 网 络 法 能 更 好 地 反映 电离 层 的 变化 ， 主 要 表现 在 正午 和 午夜 ， 而 夜间 神经 网 络 法 的 TEC 


较 低 也 与 图 5 的 结果 相符 合 。 


| 
一 40 上 
D 

Q 30 
æ 
O 20F A 
TRONS i EA E | 
0 


October 2016 (UT) 


图 8 不 同方 法 反 演 的 北纬 40.52° 西 经 3.09° 的 TEC 
Fig. 8 TEC at 40.52°N and 3.09°W derived by different methods 


我 们 提出 了 一 种 基于 人 工 神经 网 络 的 准 实时 单 站 电离 层 反 演 方法 。 在 设计 的 神经 网 络 中 ， 输 入 层 包含 


作为 初 值 的 前 一 天 的 硬件 偏差 以 及 穿刺 点 的 余 纬 、 日 地 坐标 系 下 的 经 度 以 及 高 度 角 ， 输 出 层 为 球 谐 函 数 系 
数 以 及 人 硬件 偏差 。 通 过 不 断 调整 的 权 值 ， 输 出 层 可 以 利用 梯度 下 降 算 法 根据 残 差 进行 迭代 。 因 此 硬件 偏差 
可 以 随 着 观测 值 动态 调整 并 且 TEC 也 可 以 被 准 实时 地 反 演 。 


在 一 个 实验 中 ， 位 于 北纬 30.28°? 东 经 109.46° 


的 音 


站 在 2014 年 春分 、 夏 至 、 秋 分 和 冬至 的 磁 静 日 GNSS 数据 被 用 于 分 别 采用 提出 的 方法 与 经 典 的 最 小 


二 乘 球 谐 函 数 法 进行 电离 层 反 演 。 位 于 东经 110° 子 午 线 的 18 个 站 被 作为 参考 站 通过 网 格 法 获得 硬件 偏差 


和 TEC 参考 值 。 结 果 显 示 使 用 提 昌 


的 方法 估计 的 硬件 偏差 在 整体 上 更 接近 于 参考 值 。 反 演 的 TEC 相对 于 


DE 
\\ 


最 小 二 乘 球 谐 函 数 法 更 接近 于 参考 值 。 在 另 一 个 实验 中 ，2016 年 10 月 13 日 电离 层 暴 事件 的 TEC 也 分 别 
通过 位 于 北纬 40.52° 西 经 3.09。 的 单 站 使 用 提出 的 方法 和 最 小 二 乘 球 谐 函 数 法 被 反 演 。 两 种 方法 反 演 的 TEC 


也 具有 很 好 的 一 致 性 ， 基 于 神经 网 络 的 硬件 偏差 估计 与 磁 静 日 的 硬件 偏差 估计 更 接近 。 这 些 结果 说 明了 提 


出 的 神经 网 络 方法 相 比 最 小 二 乘 球 谐 函 数 法 具有 较 高 的 精度 。 神 经 网 络 同样 适用 于 准 实时 多 站 电离 层 反 演 ， 


值得 进一步 推广 。 
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Near-real-time ionospheric TEC derivation from single station with neural 


network 


LU Weijun!:?,, MA Guanyi!:? 
(1. National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100101, China; 2. University of Chinese 
Academy of Sciences, Beijing, 100049, China) 


Abstract: In derivation of ionospheric total electron content (TEC) by global navigation satellite system 
(GNSS), compared with a multi-station observation network, single-station derivation is a flexible and convenient 
method. Based on the artificial neural network (ANN), we propose a near-real-time method to derive ionospheric 
TEC with a single station. In this method, the instrumental biases of previous period are taken as initial values and 
are adjusted with the observation data. Meanwhile, the ionospheric TEC is derived in near-real-time. In order to 
have a detailed assessment on this method, through a single station in China, the TEC during four magnetically 
quiet days is derived by the proposed method and classic least square method (LSM) with spherical harmonics, 
respectively. The references of instrumental biases and TEC are obtained by the nearby multi-station network. In 
another test, through a single station in Europe, the TEC during an ionospheric storm and the quiet days before 
and after the event are also derived by the above methods. At the magnetically quiet days, the estimated 
instrumental biases are closer to the references than LSM with spherical harmonics on the whole, and the derived 
TEC is also closer to the references. During the ionospheric storm, the TEC derived by the two methods also have 
good consistency, and the estimated instrumental biases during the ionospheric storm are closer to those at 
magnetically quiet days. The results show that the proposed neural network method has higher accuracy than LSM 
with spherical harmonics. 


Keywords: Ionosphere; Artificial neural network; Total electron content; Instrumental bias; Ionospheric storm 


